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Abstract 

By analytically constructing the matrix elements of an electron-phonon interaction for the D 
band in the Raman spectra of armchair graphene nanoribbons, we show that pseudospin and 
momentum conservation result in (i) a D band consisting of two components, (ii) a D band Raman 
intensity that is enhanced only when the polarizations of the incident and scattered light are 
parallel to the armchair edge, and (iii) the D band softening/hardening behavior caused by the 
Kohn anomaly effect is correlated with that of the G band. Several experiments are mentioned that 
are relevant to these results. It is also suggested that pseudospin is independent of the boundary 
condition for the phonon mode, while momentum conservation depends on it. 
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I. INTRODUCTION 



Raman spectroscopy has been widely used for the characterization of carbon-based ma- 
terials such as graphite, graphene, nanotubes, and nanoribbons.- Of the several bands that 
appear in the Raman spectrum, the D band at ~1350 cm -1 is a mark of the presence of 
lattice defects. 2 It is known that the formation of an electron standing wave due to inter- 
valley scattering at defects is the key to enhancing D band intensity. Recently, the D band 
has been examined intensively in relation to studies of the graphene edge.- - - In terms of 
symmetry, graphene edges can be divided into two crystal orientation categories: armchair 
and zigzag edges, and only an armchair edge causes intervalley scattering. Therefore, the D 
band is of prime importance when characterizing graphene edges. In this paper, we focus 
on graphene nanoribbons with armchair edges (armchair nanoribbons [ANRs]) to clarify the 
essential features of the D band. Several predictions that we have derived from pseudospin 
and momentum conservation for the matrix element of the electron-phonon interaction of 
the D band are useful for obtaining detailed information about the electronic and phononic 
properties of armchair nanoribbons. 

This paper is organized as follows. In Sec. [Ill we review the electronic and phononic 
properties of ANRs to provide essential background about the D band. The materials 
described in Sec. [Til are needed for calculating electron-phonon matrix elements analytically. 
In Sec. IIII| we point out several features of the D band in armchair nanoribbons. Sections IIVI 
and [V] contain our discussion and conclusion. 



II. D BAND AS FIRST ORDER PROCESS 



The electronic energy dispersion relation of ANRs is identical to that of graphene (without 
an edge).— Let k = (k x ,k y ) be the wavevector of an electron, the energy dispersion is 
given byii 

^7^1 + 4 cos 2 J^j + 4 cos J^j cos J^j , (1) 

where s (= ±1) is the band index, 7 (= 3 eV) is the hopping integral, the component k x 
(k y ) of the wavevector is perpendicular (parallel) to the edge, and a (b) denotes the periodic 
length perpendicular (parallel) to the edge [Fig. Ha)]. We use units in which a/2 = 1 and 
6/2 = 1, so that the wavevector and coordinate are dimensionless quantities. 
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FIG. 1: (a) The structure of an ANR. Carbon atoms are divided into A (•) and B (o) atoms, a is a 
lattice constant [a = 2.46 A] and b = \/3a. (b) Top: BZ of ANR. The point at (k x , k y ) = (2tt/3, 0) 
is the Dirac point. Bottom: BZ of graphene. The thick arrow denotes the phonon wavevector 
(lx,Qy) = (— 4tt/3, 0) contributing to the first order Raman process in the BZ of ANR, because of 
zone holding. 

The difference between an ANR and graphene appears in their Brillouin zones (BZs). 
The BZ of an ANR is given by k x G (0, ir) and k y G [— it/2, 7r/2) as shown in Fig. Hfb). The 
BZ of an ANR covers only one half of the graphene's BZ.— As a result, there is a single 
Dirac point in the BZ of the ANR, although there are two inequivalent Dirac points (known 
as K and K' points) in the BZ of graphene.— The BZ holding for ANRs can be understood 
more clearly by using the electronic standing wave in ANRs that is written as^ 

I / e -*e(k) 

^(..^-^c-^^ainM I | . (2) 



N 



where J (= 1, . . . , N) denotes the lattice site perpendicular to the edge [Fig. Ufa)], and x = J 
in units of a/2 = 1. The upper (lower) component of Eq. (jSJ) represents the amplitude at A- 
atom (B-atom), and the phase 0(k) is the polar angle denned with respect to the Dirac point: 
k y = (k x — 2tt/3) tan 0(k). Because G(k) is invariant throughout the intervalley scattering 
process at armchair edge,— 4>'L kxky (x) = —4>k x ,k y ( x )- Thus, (fi s kxjky (x) and (j) s _ kxtky (x) 
are identical, that is, the reflection taking place at armchair edge identifies (k x , k y ) with 
(—k x , k y ), and one half of the graphene's BZ needs to be excluded from the BZ of an ANR 
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to avoid a double counting. 

The x-component of the wavevector, k x , is quantized by the boundary condition that is 
0k( x ) = imposed at x = and x = N + 1 (i.e., at fictitious edge sites), as 



where n represents the subband index. Note that k n G (0, 7r) and that the spacing between 
adjacent k n is tt/(N + 1). This spacing is one half of the spacing obtained by imposing a 
periodic boundary condition on the wave function, which shows a double density of states 
in the reduced BZ. The double density of states in the reduced BZ is a direct consequence of 
the BZ holding and is the principal reason why the energy gap of an ANR is different from 
that of a (zigzag) nanotube obtained by rolling an ANR into a cylinder.— The reduction of 
the BZ for ANRs is not a matter of notation but a physical result since the energy gap is 
determined by the double density of states in the reduced BZ. Meanwhile, we assume that 
k y is a continuous variable. Hereafter we abbreviate Q(k n , k y ) and 4>k n ,k y by omitting k n and 
k y as O n and respectively. 

Due to the BZ holding for ANRs, momentum conservation is relaxed and a phonon 
with a non-zero wavevector can contribute to a first-order Raman process. For example, an 
intervalley phonon with a wavevector (q x , q y ) = (— 47r/3, 0) can be excited without disturbing 
the electronic state since the final electronic state (K 7 ) is identical to the initial state (K), 
as shown by the dashed arrow in the lower inset to Fig. Hfb). In graphene without an 
edge, it is considered that for a photo-excited electron at {k x , ky), only the phonon at the Y 
point (q = 0) contributes to the first order Raman process due to momentum conservation.- 
However, in ANRs, a phonon with the wavevector (q x , q y ) = (—2k x , 0) [mod tt] contributes 
to the first order Raman process.— As we will show later with the more mathematically 
rigorous method, the D band is represented as a first order, single resonance process in the 
reduced BZ, whereas it is commonly recognized as a second order, double resonance process 
in the BZ of graphene.— 1 ^ 

Vibrational displacement vectors for phonons with q = (q, 0) in ANRs can be explicitly 
constructed. First, we define the acoustic and optical branches 



k, 



(n=l,...,N) 



(3) 



n 



N + 1 
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as the basic vibrational modes in ANRs. In Eq. (jj]), the upper (lower) component represents 
the vibrational direction of an A-atom (B-atom), where (i = x,y) is the unit vector 
parallel to the i-axis. The function sin(gx) or cos(gx) represents the vibrational amplitude 
of the atom at x.— The amplitude may be normalized by multiplying each mode with the 
proper amplitude u. For the optical modes, u is nearly constant at about O.OlA, independent 
of the value of q. For the acoustic mode, u depends sensitively on the energy, which is a 
function of q and temperature. We generally have different amplitudes for acoustic and 
optical modes. However, at the point q = 2n/3 the energies of the acoustic modes become 
comparable to that of the optical modes, so the amplitudes of the acoustic and optical modes 
are approximately the same.— Then normal modes would be given by the sum of the acoustic 
and optical modes. From the basic vibrational modes given in Eq. (J3J), we construct new 
modes as 



u q {x) = sin(ga;) | + cos(gx) I | . (5) 

This mode is a superposition of the acoustic mode along the x-axis (LA) and the optical mode 
along the y-axis (TO). The mode u q (x) with q = 27r/3 corresponds to the Kekule distortion 
shown in Fig. |2]^a). In Fig. [2]^a), it is clear that the Kekule distortion does not change the 
bond angle, and that the mode is composed only of a bond stretching motion.- The mode 
u q (x) with q = can be used to represent the T point optical mode shown in Fig. Et^b) 
because the acoustic component in Eq. ([5]) disappears. We define another combination of 
acoustic and optical branches as 

v q (x) = cos(gx) I ] — sin(gx) I v ] . (6) 
\-e x J \e y J 

This mode is a superposition of the optical mode along the x-axis (LO) and the acoustic 
mode along the y-axis (TA). The mode v q (x) is given by a —90° rotation of the mode u q (x): 
v q (x) = R(—7c/2)u q (x). Thus, the mode v q (x) with q = 27r/3 is realized simply by a change 
in the bond angle, as shown in[2]^c), and the mode v q (x) with q = is the T point optical 
mode shown in Fig. |2Jd) . 

According to Tuinstra and Koenig,- the D band originates from the Kekule distortion, 
while the G band is composed of the two optical components Uo and v . The vibrational 
energy of the Kekule distortion (u 2 = 3K/M where M is the mass of the carbon atoms 1 -) is 




FIG. 2: The displacement vectors of four phonon modes, (a) Kekule distortion — it 27r /3(x), (b) 
—wo (or G y ) representing an optical phonon at the T point, (c) v 27T /s(x), and (d) vq (or G x ). 

highest among the in-plane phonons at the Dirac point, while that for v 27T /3 (uo 2 = 18H/M) 
is the lowest, because the bond stretching force constant (K) is much larger than the force 
constant for the deformation of the angle (H).— The Uq and vq energies are given by u 2 = 
3K/M + 18H/M, showing that the G band originates from both bond stretching and bond 
angle change. 2 Hereafter, we refer to the component of the G band that originates from no 
(v ) as G y {G x ) band for clarity. 

III. ELECTRON-PHONON INTERACTION 

First, we consider the electron-phonon interaction for the vector u q (x). The bond stretch- 
ing induced by a displacement vector gives rise to two distinct effects: a change in the 
nearest-neighbor hopping integral and the generation of a (deformation) potential at each 
atom. The former is known as an off-site interaction and the latter as an on-site interac- 
tion. Therefore, the electron-phonon interaction consists of off-site (i/ ff) and on-site (H on ) 
components. The matrix elements for each component are given by (see Appendix |B] for 
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derivation) 

(cj) s m \H oS (u q )\(j) s n ) 



N 



N 



sin(/c m J) cos(g J) sin(/c n J) 



j=i 

N 



2 +2sin( - + cos f- 
2 6/ V2 



— sin(/c m J) sin(gj) cos(£; n J) ) 4 sin (| + ^ sin (| j sin(fc n ) }>( 



j=i 



(7) 



and 



iV 



^|# OI1 (u g )|0*) = ( — ^sin(A; m J) cos(gJ) sin(/c n J) 



j=i 



4 sin - cos — — , 
2/ V2 6/J 



<^o> 



s s 
mn ' 



(8) 



where the matrix element of the Pauli matrix (i = x, y, z, 0) is defined as 

^i)Z=i(^ s')ai\ | , (9) 



o"o is the 2x2 identity matrix and 






(10) 



The right-hand side of the off-site component Eq. (JTj) is proportional to the pseudospin 
{<y x )mm while the on-site component depends on {<Jo)mn- The summations over J in Eq. (JTj) 
and Eq. (JB]), i.e. {2/N^2, N J=l . . .), do not vanish only when g are specific values that satisfy 
momentum conservation. So, the matrix element is determined by two factors: pseudospin 
and momentum conservation. 



Let us focus on the first order Raman process, 



(i.e., k m = k n and s' = s), that is, 



the diagonal matrix element in Eq. (J?]). By executing the summation over J for momentum 
conservation, we see that the matrix element can be non-zero if the phonon's wavevector 
q satisfies (i) q = or (ii) q = 2n — 2k n . The case (i) corresponds to the T point optical 
phonon, which is relevant to G y , while the case (ii) is relevant to the D band since q ~ 2tx/2> 
is satisfied for the electronic states near the Dirac point (k n ~ 2tt/3). 
For the D band, by setting q ~ 2tt/3 in Eq. ([7j), we obtain 



(4>n\Hos{u q ^2n/3)\(t>n) ~ H<?x) 



ss 
nn' 



(11) 
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From Eq. ([8]), the on-site component almost vanishes (0^|if on (1*^271-/3) 1 0^) — 0, so that 
only the off-site component contributes to the matrix element. The off-site component is 
proportional to {cr x ) s ^ n , which is scos0 n from Eq. flU]). The transition probability ~ cos 2 9 n 
takes its maximum (minimum) value when Q n = or tt (O n = tt/2 or — tt/2). Thus, 
the electrons at Q n = or tt contribute selectively to the D band Raman intensity, while 
the phonon emission from the electrons at B n = tt/2 or —tt/2 is negligible [Fig. |3^a)]. 
From the wave function Eq. (T5]), the photo-excited electron at 6 = (tt) corresponds to 
the bonding (antibonding) orbital. This result suggests that the bonding and antibonding 
orbitals couple strongly to the bond stretching motion, which can be also understood by 
diagonalizing the 2x2 matrix H = ^a x for a two-atom system. The eigenstates of H are 
bonding or antibonding orbitals. Since the bond stretching motion induces a change in 7 as 
7 + ^7, the off-site electron-phonon interaction is given by H Q s = 5ja x , whose expectation 
values become its maximum for the bonding or antibonding orbitals. Since the D band 
consists only of the bond stretching motion, the electrons at G n = and tt couple strongly 
to the D band phonon. 

For the G y band, by setting q = in Eq. ([7j), we obtain 

(««o)l€ = -3W (12) 

Because there is no contribution from the on-site component [(0^|if on (n o )|0^) = from 
Eq. ([!])], the transition amplitude is determined by the off-site component. It is easy to see 
that only the electrons at G n = or tt contribute to the G y band intensity. 

Comparing Eq. ffTTj) with Eq. ffT2l) . we see that the electron-phonon matrix elements for 
the D and G y bands are the same (apart from the unimportant sign change), although 
the phonon's wave vectors are distinct. A direct consequence of this fact is that if the D 
band is activated, the G y band should be activated in the same manner. For example, 
the D band intensity at armchair edge becomes comparable to the G band intensity. In 
addition, due to the fact that both the electron-phonon matrix elements are proportional to 
{cr x )n n i the D and G y band Raman intensities have the same polarization dependence: the 
Raman intensities of these phonon modes should be enhanced when the polarization of the 
incident light is parallel to the armchair edge. This polarization dependence of the D band 
has been confirmed in many experiments.- - - Several experimental groups have observed a 
correlation between the D and G bands,- - - 1 ^^ suggesting that the G x Raman intensity is 
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FIG. 3: Energy contour circle near the Dirac point, (a) The electrons at the two shaded regions 
on the energy contour circle (0 ~ or ir), have the strongest electron-phonon matrix element for 
emitting the D and G y bands. These electrons can be selectively excited by the light polarization 
setup: "in" ("sc") denotes the polarization of incident (scattered) light, (b) The electrons in the 
two shaded regions (0 ~ 7r/2 or —ir/2), have the strongest electron-phonon matrix element for 
emitting the G x band. 



suppressed at the armchair edge. Another example is the correlation between the frequency 
hardening/softening behavior of the D band due to doping (Kohn anomaly effect) and that 
of the G y band. These features of the D and G bands are further explored in the following 
subsections. 

Next, we examine the electron-phonon interaction caused by the vector v q (x) denned in 
Eq. ([6]). The corresponding electron-phonon matrix element is given by 



N 



{<j>i\H oS {v q )\<j> s n ) 



j=i 



2 sin - cos — — 
1 2/ V2 6 



— sm(k m J) sin(qj) sm(k n J) 

2 N ( 

— ^2 sin(k m J) cos(gJ) cos(k n J) J 4 cos (|j cos ^| + sin(A; r , 

(13) 



and 



N 



H n(v q )\(j) s n ) = [ — sin(fc m J) sin(gj) sin(A;„ J) 



j=i 



4 sin - cos — — 
1 2/ V2 6 



s s 
mn' 



(14) 



A notable difference between the matrix elements for u q (x) and v q {x) is that the off-site 
components are proportional to a x and a y , respectively. Moreover, the on-site interaction 
induced by u q (x) is symmetric about two sublattices, while that induced by v q (x) is anti- 
symmetric, as is evident from the matrices o"o and a z in Eq. (JSJ) and Eq. ( | 14 p . 

From Eq. (1141) the on-site component is strongly suppressed for the phonon satisfying 
q ~ or q ~ 2tt/3, and it is negligible. For q ~ 27r/3, the off-site component is also negligible 
as (4> s m \H oS (v qC:t 2iT/3)\ ( t ) n) — because cos(g/2 + n/6) ~ holds. Thus, the mode v q ~2K/z{ x ) 
does not appear in the Raman spectrum. Since 1^271-/3 W consists only of a change in 
bond angle and does not contain any bond stretching motion, the vanishing matrix elements 
are reasonably understood to show that only the bond stretching motion gives rise to the 
electron-phonon interaction. For q ~ 0, on the other hand, the off-site component can be 
non-zero as 

(^ItfoffK^oM) * -- ( 1 —~ + ) 3i(a y )i s nl (15) 

7r \m — n + p m — n — pj 

where m — n±p (= An) is an odd number and integer p is defined through q = pir/(N + 1). 
When m — n±p is an even number, the off-site element vanishes. The factor (l/nAn) on the 
right-hand side of Eq. (fl5l) . originates from momentum conservation, and a small change in 
the wave number (m — n 7^ 0) is needed to excite this G x band. A large change in the wave 
number is not important because the matrix element is suppressed by the inverse of An. 
By setting s' = s and m ~ n in Eq. (J9j), we get (cr y )f^ n ~ ssin9 n . The transition amplitude 
takes its maximum value when n = tc/2 or —71/ 2. Thus, the electrons at n = tt/2 or 
— 7r/2 contribute selectively to the G x band Raman intensity, while the phonon emission 
from the electrons at n = or 7r is negligible [Fig. [3]^b)]. In the following sections, we point 
out several consequences that can be derived from the electron-phonon matrix elements. 

A. D band splitting 

The fact that there are the two electronic states with bonding or antibonding orbitals 
(0 = or 7r) from which the D band phonon is emitted results in the splitting of the D 
band if we take account of the trigonal warping effect. The splitting width increases with 
increasing incident laser energy as 

Au D = 25 [cm- 1 ], (16) 
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FIG. 4: (a) Two principal wavevectors of the phonons (qo and q n ) contributing to the D band, (b) 
The phonon dispersion relation near the Dirac point. The energy difference between u}£>(qo) and 
u!d(Qtt) appears as two peaks in the D band. 

Since the probability that the photo-excited electron at G emits the Kekule mode is 
proportional to cos 2 [Eq. (fTTj) ]. the photo-excited electrons satisfying O ~ or ~ n, 
indicated by the shaded regions shown in Fig. [3]^a), contribute efficiently to the D band 
intensity. Thus, two phonon modes that originate from the two regions in the photo-excited 
electron contribute to the D band. We denote k x at = (0 = n) as k (k^) [see Fig. Hta)], 
where k and k n are given by solving the equation E L /2 = j\l + 2cos(k x )\ as 



2tt V3 (El\ v/3 (El 



2 



2tt . (El\ _ VS (E± 1 



3 K --2\frj) T\3^ ] 



(17) 



The ± signs in front of v / 3/8(-Bl/37) 2 represent the trigonal warping effect meaning that 
the energy contour circle is anisotropic as shown in Figs. |3]and|Ha). The corresponding 
phonon's wavevector is given by qo/ n = 2n—2k / n [Fig. H^a)]. Then we find that the difference 
between the phonon's wavevectors \q — 27r/3| (relative to the Dirac point) and \q w — 27r/3| 
is well approximated by 



For simplicity, let us assume that the phonon dispersion relation for the D band is isotropic 
about the Dirac point: uj d is a function of \q — 2n/3\ asa^dg — 27r/3|). From Fig. @Jb) it 
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is clear that the two phonon modes with go and q n have different phonon energies, which 
results in the double peak structure of the D band. 

The difference between the energy for the phonon mode with q and that for q n is ap- 
proximated by 

dun , . 

Au ° = ~&f Ag ' (19) 

where dun/dq is the slope of the phonon energy dispersion. The value of the slope can be 
estimated from the dispersive behavior of the D band: the D band frequency increases as 
the laser energy El increases. For example, Matthews et al— reported that dun/dEi for 
graphite is about 50 cm _1 /eV. A similar value has been obtained for single layer graphene.-^ 
We interpret the dispersive behavior that occurs as a result of the dispersion relation of the 
phonon mode.— Then we have 

dujp dq du D 

dE L ~ 8E L 8q ' 1 j 

By combining El/2 = j\l+2 cos(fc)| and q = 2n—2k, we have \dEi,/dq\ ~ v^T for the states 
near the Dirac point. Putting this result in Eq. ( 120]) . we have duo/dq = \f?>^{duL)/ dE L )- 
Then we obtain Eq. ffT6]) by combining Eq. f JT8|) and Eq. f fT9|) . Actual splitting can be smaller 
(or even larger) than Eq. (Tl6l) due to several factors, such as (i) the phonon dispersion 
relation for the D band not being exactly isotropic around the Dirac point, (ii) is not 
exactly limited by or n. The conditions for observing the splitting are discussed later. 



B. Light polarization dependence 

Pseudospin is the key to the characteristic polarization dependence of the D and G 
band intensities. Below we show that the Raman intensities are enhanced only when the 
polarizations of the incident and scattered light are both parallel to the armchair edge. That 
is, the D band intensity is well approximated by 

/ D (e in , e sc ) ~ icos 2 (e in ) cos 2 (e sc ), (21) 

where Oi n (0 SC ) the angle between the orientation of the armchair edge and the inci- 
dent (scattered) light polarization.— Without a polarizer for the scattered light, we obtain 
-^D(@in) - /7rcos 2 (e in ) by integrating J D (@in, @ sc ) over 6 SC . 
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To discuss light polarization dependence, let us review the interband optical ma- 
trix elements obtained in a previous paper.— The electron-light interaction is written as 
H(e) = —ev ■ A, where — e is the electron charge, v = (v x ,v y ) is the velocity operator, and 
A = Aee~ luJt is a spatially uniform vector potential. Here, hu corresponds to the incident 
light energy and e denotes the polarization. The calculated matrix elements are given in 
units of — ev^A by 

{0 m — n e even, 

--^= m-ne odd, 
7r m — n 

(CI%)|^=iC(* (23) 

The Kronecker delta 5 mn in Eq. ( 1231) shows that the y-polarized light (e y : parallel to the 
armchair edge) results in a direct interband transition (m = n). The transition amplitude 
depends on the diagonal matrix element of the pseudospin (cr x )nnr As we nave seen, Kcr^)^ 
takes its maximum value for n = or ir. On the other hand, |(o"x)^il vanishes for B n = 
±7r/2. Therefore, only the electrons at n = or ir are selectively excited by the |/-polarized 
light. From Eq. (I2"2i) . the x-polarized light (e x : perpendicular to the edge) results in an 
indirect interband transition (m — n = ±1, ±3, . . .), and the optical transition probability is 
reduced by the factor of (2/n) 2 (m — n)~ 2 , independent of the pseudospin {p y ) c ^ n . 

First, we consider a parallel configuration: both the incident and scattered lights are 
^/-polarized [see Fig. El^a)]. Since the absorption and emission of y-polarized light occur 
through direct transitions (f) v n <p c nl the phonon mode with q = or q = 2tt — 2k n fulfils the 
momentum conservation of the first order Raman process. Thus, the D and G y bands can 
be detected by the parallel configuration. The G x band is suppressed because its activation 
requires a non-zero shift in the electron wave number, and the momentum deficit cannot 
be compensated by direct optical transitions. More importantly, pseudospin enhances (sup- 
presses) the Raman intensities of the D and G y bands (G x band) because the probability 
amplitude of exciting the D and G y bands (G x band) is proportional to cos Q n (sin 0„). The 
function cos0 n agrees with the optical matrix element that is also proportional to cosB„, 
while the function sin B n is out-of-phase. The probability that the electron at B n completes 
the Raman process by emitting the D or G y band (G x band) is proportional to cos 6 0„ 
(cos 4 n sin 2 n ). We estimate the reduction in the intensity of the G x band caused by the 
disagreement of the pseudospin to be 1/5 by using the ratio of § cos 4 n sin 2 Q n dQ n (= 7r/8) 
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to § cos 6 Q n dQ n (= 57r/8). The G x band intensity is further reduced by the momentum 
deficit. 

Next, we examine a crossed configuration: the incident light is y-polarized, while the 
scattered light is x-polarized. Since x(y)-polarized light results in an indirect (direct) optical 
transition, momentum mismatch takes place: the electron cannot return to its original state 
in the valence band unless a phonon with a specific momentum is involved in the Raman 
process. The intensity is suppressed by a factor of {2/tt) 2 because of the reduction factor in 
the optical matrix element for the x-polarized light. Furthermore, momentum conservation 
for the phonon mode reduces the emission probability by the factor of (2/tt) 2 for the G x band 
[or (1/2) 2 for the D and G y bands], so that the intensity is suppressed by a factor of (2/7r) 4 
(~ 0.16). Moreover, the pseudospin for incident light does not match that for scattered 
light, and the Raman intensity is suppressed by a factor of 1/5 in the crossed configuration. 
Thus, the D and G y bands intensities are suppressed by a few percent compared with the 
intensities in the parallel configuration. 

Finally, we study a perpendicular configuration: both the incident and scattered lights 
are x-polarized [see Fig. EJ^b)]. The Raman intensities of the D, G y , and G x bands for this 
configuration are suppressed by (2/7r) 4 at least, due to the reduction factor for the optical 
matrix elements. The pseudospin for the x-polarized light agrees only with that for the 
G x band. Thus, the D and G y bands are not detectable in the perpendicular configuration 
because of the further reduction by a factor of 1/5. 

When the polarization of the incident light is perpendicular to the armchair edge, it 
is common for the D band to show a residual intensity.-^"- The ratio of the minimum 
D band intensity (observed when the incident light polarization is perpendicular to the 
edge) to the maximum intensity (observed when the incident light polarization is parallel 
to the edge) is 0.12 ~ 0.2 depending on the experiments. Several authors have suggested 
that the presence of the residual D band intensity is correlated with the irregularities that 
consist of armchair segments having angles of ±60° with respect to the armchair edge.-^ 2 ^ 
Because their arguments assume that the irregularities contribute independently to the total 
D band intensity as ^(©in) — acos 2 (6 in ) + 6{cos 2 (6 in + 60°) + cos 2 (6 in — 60°)}, there is 
a possibility that this assumption is not plausible when the electron wave function extends 
along the armchair edge. We provide a different interpretation of the residual D band 
intensity. When the momentum conservation that we have used is not physically legitimate 
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because of the boundary condition of the phonon mode (see Appendix [B] for a relevant 
discussion), we can only rely on the pseudospin. Then we may expect the residual D band 
intensity to be about 0.2 (1/5). Even for an armchair edge without irregularity, the D band 
exhibits a residual intensity depending on the boundary condition of the phonon mode. 

Without a polarizer for the scattered light, the D band polarization dependence is ex- 
pressed as 

/ D (©in) = acos 2 (9 in ) + 6sin 2 (9 in ) + c, (24) 

where a, b, and c are parameters. When both the momentum conservation and pseudospin 
are physically legitimate, b ~ 0. When we can only rely on the pseudospin, we expect 
b/a to be 1/3 by using the ratio of § cos 2 6 sin 2 QdQ (= tt/4) to ^cos 4 6d9 (= 3vr/4). 
Figure [5] shows the polar plots for these two cases. If the residual D band intensity decreases 
by putting a polarizer for the scattered light, such a signal would evidence the effect of 
the pseudospin. The 6 in independent c term may originate from a point defect or on-site 
deformation potential (see Appendix [A} and obscures the D band polarization dependence. 



0° 




FIG. 5: The polar plot for the D band intensity. The parameters for the solid curve are a = 1, 
b = 1/3, and c = 0, while those for the dashed curve are a = 1 and b = c = 0. 

C. Kohn anomaly and transport 

In ANRs, because of the conical shape of the energy spectrum, momentum conservation 
tells us that a vertical (direct) electron-hole pair can couple to the T point phonon. Fur- 
thermore, due to the BZ holding, the direct electron-hole pair can also couple to the D 
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band phonons. As a result, the G and D bands in Raman spectra can exhibit a similar 
Kohn anomaly effect. The probability amplitude for the process whereby the D and G y 
bands change into a vertical electron- hole pair is given by using |$^) = <t z |0^), a x a z = —icr y , 
Eq. dHj, and Eq. ([12]), as 

(25) 

Since {u y )n n = sin0 n , electrons at Q n = n/2 or — 7r/2 contribute to the Kohn anomaly. 
Equation (125]) clearly shows that the behavior of the Kohn anomaly for the D band is 
correlated to that for the G y band: if the D band exhibits hardening/softening due to the 
Kohn anomaly effect, the G y band also exhibits hardening/softening. Recently, Zhang and 
Li^ have observed a Kohn anomaly effect for the G band at an armchair graphene edge. 
Interestingly they also found a correlation between the 2D band and the G band. 

It is interesting to note that the pseudospin for electron-hole pair creation is identical to 
the matrix elements of the velocity operator along the armchair edge, which is written as 
\v y 1 0°) = S mn VF{a y )^ n .— Thus, the Kohn anomaly effect is fundamentally related to the 
velocity or current behavior. If the velocity is suppressed, then the Kohn anomaly effect 
is suppressed as well. The velocity perpendicular to the armchair edge v x is suppressed 
by electron reflection, while the velocity along the edge v y is not. As a result, only the D 
and G y bands can experience a strong Kohn anomaly effect.— In the presence of impurities, 
the electronic velocity (current) along the edge (ribbon) might also be suppressed due to 
the scattering. However, in metallic ANR, because Berry's phase can protect the electronic 
current along the armchair edge from decaying,— the Kohn anomaly effect for the D and G y 
bands should be robust against impurity potentials. We consider that the D and G y bands 
undergo a strong Kohn anomaly effect even when we take account of the effect of impurity 
potentials. 

IV. DISCUSSION 

Here we discuss the possible factors that obscure the double- component feature of the D 
band proposed in Sec. IIII Al First, the resonance effect suppresses the intensity of one of the 
two components. Suppose that the electron at ko is resonant with the incident laser energy. 
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Then the Raman intensity for the process in which the photo-excited electron at ko emits 
a phonon with momentum go is enhanced as compared with the Raman intensity involving 
a phonon with momentum q n . Such a resonance effect is apparent in thin ANRs, whereas 
it is negligible for thick ANRs if the ANR width (W) exceeds a critical value. The critical 
W is estimated from the condition whereby the difference between fc ~~ 27r/3 and 2tt/3 — k n 
caused by the trigonal warping effect [Eq. (IT?]) ] is larger than the spacing of each sub-band, 
tt/(N + 1). This condition leads to W > 8^/El) 2 nm: W > 13 nm when E L = 2.3 eV. 
Second, for thick ANRs satisfying W > S^/El) 2 nm or for armchair edges of graphene, the 
energy gap is smaller than the energy of the D band. As a result, the width of each peak can 
be wide enough to obscure the splitting structure since the D band undergoes a strong Kohn 
anomaly effect (Sec. IIII Cp . To suppress the broadening of each peak caused by the Kohn 
anomaly effect, the Fermi energy needs to be placed away from the Dirac point (|-Ef| > 0.1 
eV). Third, the parallel configuration is the most suitable polarization setting with which 
to observe the splitting of the D band because the probability cos 6 G exhibits prominent 
peaks at = or it. To observe the splitting of the D band, it is important to recognize 
that the results obtained in Sees. IIII Al IIII Bl and IIII CI are closely correlated. Among the 
conclusions obtained in Sees. IIII IIII B] and IIII CI we think that the light polarization 
dependence of the D band Raman intensity is the most direct conclusion derived from the 
pseudospins of the electron-phonon and optical matrix elements. Indeed, many studies have 
reported the D band polarization behavior.- - - 1 ^ However, there is a small possibility that 
only the pseudospin of the optical matrix element is the result of the observed polarization 
dependence for the D band if there is a strong resonance effect.- 5 A strong piece of evidence 
for the pseudospin of the electron-phonon matrix element is provided by the splitting of the 
D band. 

A multicomponent D band has been seen in several sets of Raman data obtained at the 
graphene edge. Ferrari et alr& and Gupta et air- reported that for graphite, the D band can 
be well fitted by a doublet with two broad components, although a single-layer graphene 
edge produces a narrow single- component D band. Our theory relates to the armchair edge 
of a single-layer graphene. However, the results that we have obtained for ANRs may be 
applicable to single-wall carbon nanotubes because defects in carbon nanotubes result in 
the formation of a standing wave. Suzuki and Hibinc 27 reported that there is a correlation 
between the behaviors of the D and G bands under doping. Moreover, they observed D 
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band splitting when the excitation laser energy was at its largest (532 nm), while such a 
splitting was not seen when they used lower excitation energies, 633 nm and 785 nm. These 
observations are consistent with our results, that is, Kohn anomaly effects for the D and G 
bands are correlated and the D band splitting increases with increasing El- We note that 
work has already been published on the splitting of the D band in carbon nanotubes,— ^ 
in which the authors propose a completely different mechanism from ours (trigonal warping 
effect for phonon dispersion). 

V. CONCLUSION 

By constructing the matrix element of the electron-phonon interaction for the D band 
analytically, we have shown that the D band in an ANR is represented as a first order, 
single resonance process in a reduced BZ, as well as the G band. This explains clearly and 
naturally the fact that the D and G bands show a similar intensity in many experiments 
for armchair edge. The matrix element is proportional to the pseudospin ((J x )nni an d the 
D band couples selectively to the two pseudospin states (u x )nn = corresponding to 
the bonding and antibonding orbitals. The pseudospin is the origin of the branches of 
phenomena associated with the D band in ANRs, as shown in Sees. IIII Al IIII Bl and IIII CI 
The polarization dependence of the D band is consistent with the previous experimental 
reports. Direct evidence for the pseudospin of the D band is provided by the D band 
splitting at the armchair edge, which has not been previously reported. 
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Appendix A: Generalized Kekule distortion 



When the amplitudes of the acoustic and optical modes are not exactly the same, we 
may consider the extension of Eq. to the general case as 

(1 - e) sin(gx) j ^ ] + (1 + e) cos(ga;) I &y . (Al) 
\e x J \-e y J 

In this Appendix, we show how the analysis is affected by the deviation e. 

The generalized displacement vector is written as u q (x) + eu- q (x), and the corresponding 

electron-phonon matrix elements are 



m I 



|#on/offK + = (0m I ^on/off («,) Wn) + e(</>ml#on/off (™- ? ) \ <j>' n ) . (A2) 

Because the last term is given by replacing q with — q in Eqs. ([7]) and (jSj), we see that 
the off-site component is suppressed for q ~ 2ir/3, e(4>^ l \H g(u_ q ^_27 L )\(j)^} ~ 0, and that 
the correction to the matrix element of Eq. ( fTTT) arises only from the on-site component, 
e (0^ | i?o n (u_ 9 ^_27r )|</>*) ~ — 3e(oo)mrr Thus, the total electron-phonon matrix element for 
the deformed Kekule distortion is given by 

3#off (°z)mn - 3e &m(o-o)mn> (A3) 

where g on (g s) denotes the coupling constant of the on-site (off-site) deformation potential. 
Since the intensity is relevant to the square of the matrix element, the correction can in- 
clude the crossing term — 18eg on g o s(a x )^ n (a )^ n and the order e 2 term 9e 2 g 2 n ((a" );^J 2 . The 
crossing term does not contribute to the intensity because it vanishes after the integral over 
B n . Thus, the leading contribution is given by the order e 2 . The effect of the deviation is 
negligible when |e| <C \g s/g n\, where \g s/g n\ ~ 1/3 was suggested by the result of density 
functional theory— Note also that the pseudospin is proportional to er , which shows that 
the deviation is not relevant to the polarization dependence of the D band intensity but 
relevant to the c parameter in Eq. (|24|) or residual D band intensity. 

Using Eq. (IA3j) . a quick comparison can be made between the predictions of our model 
on which the D band is described as coming from a Kekule distortion and other models 
where the D band is described solely as a combination of intervalley optical phonons. The 
consequence of other models may be obtained by setting e = 1 in Eq. (IA3I) . Then, the 
matrix element is dominated by the on-site component and the pseudospin a x is suppressed. 
As a result, the D band does not exhibit the polarization dependence. 
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Appendix B: Derivations of Eqs. ([7]) and ([8]), and boundary condition for phonon 



Below we give the derivation of the electron-phonon matrix elements. Note that u q (J) in 
Eq. (jSJ) consists only of the optical mode at the fictitious edge site of J = 0. This means that 
Eq. (jSJ) results from a specific choice of phonon boundary condition. To make the derivation 
general enough to cover the possible boundary conditions, we first generalize u q (J) of Eq. (J5]) 
by inserting a phase shift of <p as 

sin(gj + cp) \ x + cos(g J + (p) I By . (Bl) 
\e x J \-e y J 

For ip = 0, Eq. (jEU) reproduces Eq. (JSJ). For (p = ir/2, Eq. (jET) leads to 

u q {J) = cos(gJ) f Gx - sin(gj) f &v . (B2) 

This u q {J) consists only of the acoustic mode at J = 0, which is contrasted with that 
u q (J) of Eq. (jSJ) that consists only of the optical mode at J = 0. 36 The phase (p is not 
very meaningful for a periodic system without an edge, such as a nanotube, however, for 
a nanoribbon with an edge, phase (p should be taken into account because it is possible 
that the phonon boundary condition could be sensitive to the situation of the armchair 
edge. It turns out that tp changes the momentum conservation slightly, while it does not 
modify the pseudospin structure of the matrix element. Thus, an analysis based on the 
pseudospin seems plausible, while the results based only on the momentum conservation 
might be physically fragile. Similarly, the v q (J) mode can be generalized as 

cos(qJ + <p)l Gx I - sin(gj + (p) [ ^ I . (B3) 
\-e x J \e y J 

In addition, the physics of graphene edge is, in some sense, an understanding of the boundary 
conditions at edges. The importance of the boundary conditions for electrons has been widely 
recognized, while those for the phonon mode and for electron-phonon interactions are yet 
to be understood. 

The off-site electron-phonon matrix elements for the displacement vector u q (J) of 
Eq. (1B1I) are written as 

N , 

(<|F off KM> = £ {5h + jG + ^ J+1 + Shj^j + ShjG-^) . (B4) 

j=i 
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On the right-hand side, we have defined the translational operators, 

/ e ikyb Q \ fi \ 

' („ :) ' ( ' (B5) 

and the 2x2 matrices of the deformation potential as 

Shj = 




6hj = 
where 

5 7 i(J)=[<(J)-<(J)]-e 1 , 

8j 2 (J) = [<( J - 1) - <( J)] -e 2 , (B7) 

^(J)= [«f(J+l)-uf(J)] .e 3l 

with Ug(J) = sin(gj + y2)e x + cos(gJ + v?) e ?/ and u^(J) = sin(gj + y^e^, — cos(gJ + v?) e ?/ 
from IBll In Eq. ( 1B7I) e a (a = 1,2,3) are the dimensionless unit vectors pointing from an 
A-atom to the nearest-neighbor B-atoms: 

ei = e y , e 2 = — —e x - -e y , e 3 = — e x - -e y . (B8) 



By inserting Eq. f ]Blj) into Eq. ( 1B7j) . we rewrite Eq. (1B6|) as 

= — 2<7a; cos(gJ + v?), 



2cr x cos J + if + |j ■ e 3 , (B9) 
5hj = 2a x cos (q J + <p - | ) e y ■ e 3 , 



where we have defined 



e q = sin ( I ) e a . - cos (| ) e y . (BIO) 



It is easy to show that e q satisfies 



' q 7f s 

e q ■ e 3 = sin (- + -)= e_ g • e 2 . (Bll) 
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Finally, by putting Eq. (1B9j) into Eq. f ]B4j) . we obtain 



j=i 



(CI#off(« 9 M) 

AT t 

= H [^m,j] { 2 Te 9 ■ e 3 cos (qj + ip + |j <^ J+1 - 2a, cos (gj + tp) <t> s n>J + 2T- 1 e g • e 3 cos (q J + <p - |j 

.7=1 
/V . 

[0m,J {-2cos(gJ + ip) a x + K s n e q • e 3 cos (|) 4> s n J ~2sm{qJ + tp)e q ■ e 3 sin (|) [T</4, 7+1 - T~ 1 <f> s n j_ 
2 N \ 

— sin(/c m J) cos(q J + </?) sin(/c„ J) J 2 + 2e q ■ e 3 cos (J^J 

^ sin(/c m J) sin(<?J + 99) cos(k n J) J 4e, • e 3 sin (|) sin(fc„) I {<J x )m n i 



(B12) 



.7=1 



where the matrix T (T 1 ) is defined by T = o^G^ (T 1 = cr z G ). Note that the energy 
eigen equation for the electron in an ANR is written in terms of these T, T _1 , and K* 



Tte J+1 + Ktf-j + T 



J n,J-l 



0. 



(B13) 



The energy eigen equation has been used to obtain the third line in Eq. flB12p . Note also 
that the small term, e s kn fe /7, has been omitted to obtain the last line. By using Eq. (T5]), we 
obtain the following equation, 



\ s 



(B14) 



which has been used to get the last line. We reproduce Eq. ([7]) by putting ip = in Eq. (IB12j) . 
When ip = 7r/2, the momentum conservation in Eq. (1B12I) changes slightly. It is important 
to mention that the pseudospin is independent of the choice of the value of <p. Because the 
matrix element for the generalized v q (J) can be obtained by replacing u q (J) with v q (J) in 
the above calculation, we omit the derivation of the matrix elements for v q (J). 
The on-site deformation potential is written as 



Sh, 



Shj 
5h% 



(B15) 



where 



Sh 1 , 



u B q {J) • ei + uJ(J - 1) • e 2 + u*{J + 1) • e 3 , 
- «( J) • e x + uf(J - 1) • e 3 + <( J + 1) • e 2 } 



(B16) 
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The matrix element is obtained by inserting the above 5hj into Eq. (1B4j) and setting 8h^ = 0. 
By putting Eq. (jET) into Eq. (jEl6j) . we obtain for J = 2, . . . , N - 1 that 

= 4a cos(gJ + <p) sin | cos + |) > ( B17 ) 

which reproduces Eq. ([8]) when tp = 0. The calculation of the on-site deformation potential 
at the boundary, 5h\ and 5hw, requires careful treatment because there is no site at J = 
and J = N + 1. In general, we have a boundary deformation potential, and we cannot use 
Eq. (IB17P to estimate the boundary deformation potential. This fact can be understood in 
Fig. MJd) by looking at the relative displacement of the A (B) atom at the boundary site with 
J = 1 in relation to that of the nearest-neighbor B (A) atom at J = 2. For example, we have 
Shi(u q ) = diag(l/2, 1/2) even if q = 0. This is not expected from a simple application of 
Eq. ( IB 171) to the boundary site at J = 1. The effect of the boundary deformation potential 
on the electronic state may prove important when we consider localized states. 
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as the 2D band (overtone of D band). 

Note that the BZ for a phonon is given by one half of the graphene's BZ, q £ [0,ir], as well as 
the ANR's BZ for an electron. 

The following analysis is valid even though we take into account the difference of the amplitudes 

of the acoustic and optical modes. See Appendix [A] for detailed discussion of this point. 

At the graphene edge, it is difficult to expect a resonance effect, such as the Van Hove singularity 

for electrons and phonons. The polarization dependence of the D band at a graphene edge 

supports the pseudospin dependence of the electron-phonon matrix element. 

The radial-breathing like mode described by Zhou and DongSi corresponds to the acoustic mode 

with the form cos(qJ)(e x ,e x ) with q = tt/(N + 1). 
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